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ABSTRACT 

A  discrete-time  discrete-space  search  model  is  con- 
sidered in  which  an  observer  employing  an  idealized  de- 
tection device  is  searching  for  a  uniformly  distributed 
stationary  target.   The  model  is  formulated  as  a  discrete- 
time  counting  process,  called  the  search  process,  which 
under  weak  additional  conditions  is  uniquely  determined 
by  a  sequence  of  probabilities.   Formulas  for  the  time- 
to-detection  and  the  detection  rate  of  a  search  are 
derived  in  terms  of  the  parameters  of  the  search  process, 
and  are  applied  to  two  special  types  of  searches,  the 
systematic  search  and  the  random  search.   Using  these 
search  types  as  boundary  cases  a  purposeful  search  is 
defined,  and  sufficient  conditions  on  the  sequence  of 
probabilities  are  established  for  the  purposeful  search. 
Possible  extensions  of  the  search  process  to  less  re- 
stricted models  are  indicated. 
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1.   INTRODUCTION 

Systematic  research  in  the  field  now  commonly  known 
as  search  theory  was  begun  during  World  War  II  when  the 
Anti-Submarine  Warfare  Operations  Research  Group  (ASWORG) 
of  the  U.  S.  Navy  studied  the  German  submarine  threat  to 
allied  shipping  and  developed  methods  and  tactics  to 
counter  this  threat.   The  major  work  done  by  this  organi- 
zation is  covered  in  Koopman's  book  Search  and  Screening 
[Ref.  1]  which,  although  now  more  than  25  years  old,  is 
still  the  fundamental  work  on  the  subject. 

In  the  years  since  then  the  field  of  search  has  grown 
and  expanded.   A  number  of  specific  problems  have  been 
formulated  and  solved,  new  areas  of  application  have  been 
added  (e.g.,  mineralogical  surveying,  commercial  fishing), 
new  methodologies  have  been  applied  (e.g.,  game  theory, 
decision  analysis),  and  researchers  from  other  disciplines, 
such  as  electrical  engineers,  have  become  involved. 
Enslow's  bibliography  with  abstracts  [Ref.  2]  and  Dobbie's 
[Ref.  3]  and  Pollock's  [Ref.  4]  articles  provide  an  over- 
view of  this  growth.   Despite  this,  "the  field  still  ap- 
pears to  be  relatively  unstructured"  and  "very  little 
general  theory  has  evolved"  as  Enslow  remarks. 

The  present  paper  approaches  search  theory  from  a 
point  of  view  which  apparently  has  not  been  taken  before. 


1Enslow  [Ref.  2],  p.  177 


In  doing  so,  the  aim  is  not  to  solve  another  specific 
problem,  i.e.,  to  find  an  optimal  search  method  for  some 
particular  context,  but  rather  to  describe  the  basic  and 
essential  structure  of  all  "meaningful"  or  "purposeful" 
searches  in  the  hope  that  this  may  be  a  way  to  get  closer 
to  a  unifying  theory  of  search. 


2.   THE  SEARCH  MODEL 

Throughout  this  paper  the  following  search  model  is 
considered: 

One  observer  is  searching  for  one  stationary 
target  which  is  equally  likely  to  be  in  any  one 
of  n  squares  which  constitute  the  search  area. 
At  each  time  step  the  observer  searches  one 
square;  if  the  target  is  in  that  square  it  is 
detected  and  the  search  ends;  if  the  target  is 
not  in  that  square  the  observer  knows  that  it  is 
not  in  that  square  and  goes  on  to  search  another 
(possibly  the  same)  square.   The  observer  is 
not  constrained  to  search  the  squares  in  any 
particular  order. 

The  model  assumes  that  the  observer  uses  a  definite 

2 
range  law  detection  device   which  has  zero  probability 

of  false  alarms  .   This  assumption  implies  that  any 

observer  who  wants  to  detect  the  target  would  never 

want  to  search  the  same  square  more  than  once.   It  is 

not  unreasonable,  however,  to  assume  that,  unwillingly 

and  unknowingly ,  he  may  search  squares  more  than  once 


2 

A  definite  range  law  detection  device,  commonly 

called  a  "cookie  cutter"  device,  detects  a  target  within 
its  range  with  probability  one. 

3 
A  false  alarm  is  the  report  of  a  detection  when  in 

fact  the  target  is  not  present. 


due,  for  instance,  to  navigational  errors,  limited  memory, 
etc.   Hence  the  actual  path  of  the  observer  through  the 
search  area  is  determined  by  chance.   Of  course,  when  the 
observer's  detection  device  is  not  of  the  definite  range 
law  type,  searching  a  square  more  than  once  may  be  a  sen- 
sible course  of  action. 

Clearly  the  model  deviates  very  much  from  the  reality 
of  traditional  naval  searches  such  as  a  raider  searching 
for  merchant  ships,  destroyers  hunting  a  submarine,  or 
SAR  (search  and  rescue)  units  trying  to  locate  the  crew 
of  a  downed  aircraft .   As  far  as  modern  methods  of  anti- 
submarine warfare  are  concerned  (e.g.  those  which  employ  ASW 
helicopters  with  sono  buoys)  it  may  be  a  reasonable  first 
approximation;  and  there  are  situations  which  it  describes 
rather  accurately;  consider,  for  instance,  a  student  trying 
to  find  a  formula  which  he  knows  is  written  on  one  of  the 
many  pieces  of  paper  scattered  across  his  desk. 


3.   MATHEMATICAL  FORMULATION  OF  THE  SEARCH  MODEL 

The  model  under  consideration  can  be  described  math- 
ematically in  a  number  of  ways.   Since  the  position  of 
the  target  and  in  most  cases  also  the  path  of  the  observer 
are  chance  dependent,  any  description  will  necessarily 
be  probabilistic  (rather  than  analytic)  in  nature.   A 
discrete-time  counting  process  has  been  chosen  here 
because  it  is  felt  that  this  particular  formulation  has 
intuitive  appeal  and  exposes  the  essential  features  of 
a  search  quite  well.   It  will  be  developed  in  three 
stages. 

3.1  THE  AREA  ACQUISITION  PROCESS 

Suppose  for  the  moment  that,  unknown  to  the  observer, 
the  target  is  not  present  in  the  search  area,  and  assume 
that  the  search  begins  with  time  step  one.   Define  A(0)  e  0, 
and  for  m  =  1,2...   let  A(m)  be  the  number  of 
different  squares  which  the  observer  has  searched  by 
time  step  m.   Since  at  each  time  step  the  observer 
searches  exactly  one  square,  one  square  is  searched  by 
time  step  one  and  hence  A(l)  =  1.   At  each  successive 
time  step  the  square  searched  may  be  either  one  which  was 
not  searched  previously  (a  "new"  square)  or  else  one 
which  was  searched  previously  (a  "used"  square),  until 
all  n  squares  have  been  searched.   Then  only  used  squares 
can  be  searched.   Thus  A(m),  m  =  1,2,...,  are  random 
variables  which  take  on  the  values  1,2,...,  min{m,n} . 
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The  sequence   <A(m) ,  m  =  0,1,...  >   will  be  called 
the  area  acquisition  process  and  the  search  of  a  new 
square  an  acquisition  event .   The  k —  event  time ,  S(k), 
is  the  number  of  the  time  step  at  which  the  k —  acqui- 
sition event  occurs,  i.e. 

S(k)  =  min(m|A(m)  =  k},  k  =  l,...,n+l. 
The  inter-event  times  T(k)  are  defined  by 

T(l)  =  S(l) 

T(k)  e  S(k)-S(k-1),  k  =  2,...,n+l. 

Since  A(l)  =  1,  thus  S(l)=l  and  T(l)  =  1;  also,  since 

st 
an  (n+1) —  acquisition  event  can  never  happen  (recall 

that  there  are  only  n  squares),  S(n+1)  and  T(n+1)  are 

both  infinite. 

In  the  acquisition  process  the  fundamental  relationship 
between  the  number  of  events  and  the  event  times  in  a 
counting  process  takes  on  the  form 

{S(k)  <  m}  iff  {A(m)  >_  k}. 

To  derive  the  distributions  of  the  inter-event  times 

consider  the  random  variables  X(m)  defined  by 

'1,  if  an  acquisition  event  occurs  at 
X(m)=.    time  step  m 

"0,  otherwise 

for  all  positive  integers  m. 

For  most  practical  searches  one  would  expect  the 

probability  of  searching  a  new  square  at  time  step  m 

to  be  a  function  of  the  number  of  different  squares  searched 

before  time  step  m  (the  greater  the  number  of  squares 

searched  the  fewer  new  squares  are  left  and  the  more 

likely  an  old  square  is  searched  again) ,  but  not  a  function 


of  the  time  step  itself.   Obviously,  exceptions  are 
possible  (the  crew  of  the  observer  may  get  tired  and 
make  navigational  errors  after  some  time,  no  matter 
how  many  different  squares  have  been  searched).   Moreover, 
other  factors  may  also  influence  this  probability:  the 
weather,  the  time  of  day,  the  observer's  navigator,  etc. 
These  factors,  however,  are  neglected  here,  and  it  is 
postulated  that  the  probability  is  a  function  only  of  the 
number  of  different  squares  searched  so  far,  i.e. 

P[X(m)  =  l|A(m-l)  =  k]  =  a(k),  m  =  1,2,..., 
where  a(k)  is  a  parameter  which  depends  on  k  only. 
Clearly, 

a(0)  =  P[X(1)  =  1|A(0)  =  0]  =  1 
since  the  first  square  searched  is  always  a  new  square, 
and 

a(n)  =  P[X(m)  =  l|A(m-l)  =  n]  =  0 
because  no  new  square  can  be  searched  after  n  squares 
have  been  searched.   For  intermediate  values  of  k 

0  <_  a(k)  <_  1, 
since  the  a(k)'s  are  probabilities. 

The  distribution  of  the  inter-event  times  can  now  be 
derived;  for  any  k  =  l,...,n+l  and  any  j  =  1,2,..., 
P[T(k)  >  j]  is  the  probability  that  at  all  of  the  j  time 

steps  immediately  following  the  time  step  at  which  the 

st 
(k-1) —  event  has  occurred  no  event  will  occur.  By  the 

definition  of  a(k),  however,  there  is  probability 

[l-a(k-l)]  that  no  event  will  occur  at  each  of  the  j  time 
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steps;  hence  P[T(k)  >  j]  =  [l-a(k-l)]^.   Thus  it  follows 
that  T(l) , . . . ,T(n+l)  are  independent  and  that  T(k)  has 
the  geometric  distribution  with  success  probability 
p  =  a(k-l).   In  particular,  T(l)  =  1  and  T(n+1)  =  +«>, 
which  agrees  with  previous  results. 

3.2  THE  TARGET  VARIABLE 

The  target  and  its  detection  which  have  been  neglected 
so  far  are  now  taken  into  account  by  introducing  the 
detection  event  into  the  area  acquisition  process.   This 
event  is  bound  to  occur  simultaneously  with  an  acquisition 
event  because  the  target  can  be  detected  only  when  a  new 
square  is  searched,  and  is  equally  likely  to  occur  at  any 
one  of  the  n  acquisition  events  since  the  target  is 
uniformly  distributed  over  the  n  squares.   Stated  formally 
the  detection  event  occurs  at  the  Z —  acquisition  event 
where  the  random  variable  Z,  called  the  target  variable , 
is  independent  of  the  acquisition  process  and  has  proba- 
bility mass  function  P[Z  =  k]  =  — ,  k  =  l,...,n. 

3.3  THE  SEARCH  PROCESS 

The  area  acquisition  process  and  the  target  variable 
together  form  the  search  process .   Its  explicit  definition 
is  given  here  as  a  summary  of  the  results  of  this  section. 

The  search  process  is  a  pair  {<A(m),  m  =  0,1, ...>,Z} 
where 

(1)  the  area  acquisition  process  <A(m),  m  =  0, !,...> 


11 


counts  the  number  of  acquisition  events  in  any 
interval  (0,m];  its  inter-event  times  T(k), 
k  =  l,...,n+l,  are  such  that 

(a)  T(l) , . . . ,T(n+l)  are  independent,  and 

(b)  T(k)  has  the  geometric  distribution  with 
success  probability  p  =  a(k-l)  where 

(c)  the  parameters  a(k),  k  =  0,...,n,  satisfy 
the  relation  0  =  a(n)  <_  a(k)  <_  a(0)  =  1; 

(2)  the  target  variable  Z  with  probability  mass 
function  P[Z  =  k]  =  — ,  k  =  l,...,n,  is  indepen- 
dent of  the  acquisition  process;  and  which 

(3)  terminates  with  the  detection  at  the  time  of 
the  occurrence  of  the  Z —  acquisition  event. 

The  search  process  constitutes  the  mathematical 
formulation  of  the  proposed  search  model  under  the 
assumption  that  the  probability  of  searching  a  new  square 
depends  only  on  the  number  of  different  squares  searched 
so  far.   It  is  important  to  note,  and  obvious  from  the 
definition,  that  the  search  process  is  uniquely  speci- 
fied by  a  sequence   <ct(k),  k  =  0,...,n>  . 
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4.   TIME-TO-DETECTION  AND  DETECTION  RATE 

Two  quantities  of  particular  interest  in  a  search 

are  the  time-to-detection  and  the  detection  rate. 

The  time-to-detection,  D,  is  the  number  of  the  time 

step  at  which  the  detection  occurs;  in  terms  of  the 

search  process  it  is  the  number  of  the  time  step  at  which 

the  Z —  acquisition  event  occurs,  i.e. 

D  =  S(Z). 

Its  distribution  function,  denoted  by  FD,  can  be  expressed 

as  follows : 

FD(m)  £  P[D  <  m] 

=  P[S(Z)  <_  m] 

n 
=  I      P[S(Z)  <_   m|Z  =  k]  •  P[Z  =  k] 
k=l 

n 
=  I      P[S(k)  ±  m]  •  P[Z  =  k] 
k=l 

n 
=  l      P[A(m)  >  k]  •  P[Z  =  k] 
k=l 

1   n 
*=  i-  l      P[A(m)  >  k] 

n  k=l 
=  1-  EA(m).  (1) 

In  a  discrete-time  model  the  detection  rate ,  d(m), 
is  defined  to  be  the  probability  that  detection  occurs 
at  time  step  m  given  it  has  not  occurred  before, 
d(m)  e  P[D  =  m|D  >  m-1] , 
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4 
for  all  integers  m  such  that  P[D  >  m-1]  >  0  .   d(m) 

can  be  written  in  terms  of  F~  as 

d(m)  =  P[D  =  m|D  >  m-1] 

P[D  =  m] 
P[D  >  m-1] 

FD(m)-FD(m-l) 
l-FD(m-l) 

or,  using  equation  (1),  as 

w  ^    EA(m)-EA(m-l) 

d(m)  =  —  57  / — H^ — -  (o's 

n-EA(m-l)  \^J 

Consider  now  the  quantity  EA(m).   By  definition, 

A(m)  is  the  number  of  different  squares  searched  by  time 

step  m  which  is  the  same  as  the  number  of  time  steps  up 

to  and  including  m  at  which  a  new  square  is  searched; 

hence 

m 
A(m)  =  E   X(i)  (3) 

i  =  l 


and 


m 

EA(m)  =  E  Z   X(i) 
i-1 

m 
=  I      EX(i) 
i=l 


m 
=  i     P[x(i)  =  l],  (4) 

i=l 


since  X(i)  is  a  Bernoulli  random  variable 


This  condition  is  always  satisfied  unless  a(k)  =  1 
for  all  k  =  1, . . . ,n-l. 


1H 


P[X(i)  =  1]  can  be  expressed  using  the  definition  of 

a(k)  as 

P[X(i)  =  1]  =  l    P[X(i)=l|A(i-l)=k]  •  P[A(i-l)=k] 
k 


=  E  a(k)  •  P[A(i-l)  =  k] 
k 


=  EoCA(i-l)], 
which  when  combined  with  equation  (4)  yields 


m 
EA(m)  =  Z   Ea[A(i-l)],  m  =  1,2,...   .        (5) 
i=l 


Equation  (5)  can  now  be  used  to  write  equation  (2) 
in  a  different  way.   Since 


m  m-1 

EA(m)-EA(m-l)  =  I   Ea[A(i-l)]  -  Z   Ea[A(i-l)] 

i=l  i=l 


=  Ea[A(m-l)], 

thus 

,  s  o  Ea[A(m-l)]      m  =  1,2,...,  (6) 

(Um;    n-EA(m-l)   »    n  >  EA(m-l).  K° } 

Finally  the  distribution  of  the  time-to-detection 

can  be  expressed  in  terms  of  the  detection  rate.   From 

P[D  >  m]  =  P[D  ¥■   m,  D  >  m-1] 

=  P[D  j   m|D  >  m-1]  •  P[D  >  m-1] 

=    Q-P[D   =  m|D    >    m-1]}    •    P[D    >   m-1] 

=    [l-d(m)]  •    P[D   >    m-1] 

=    [l-d(m)]  •    [l-d(m-l)]    •    P[D   >   m-2] 
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P[D  >  m]  =  [l-d(m)]  •  ...  •  [l-d(l)]  •  P[D  >  0] 

m 
=  n  [l-d(i)] 
i=l 


it  follows  that 


m 


FD(m)  =  1  -  n   [l-d(i)]   .  (7) 

i  =  l 
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5.   SYSTEMATIC  SEARCH  AND  RANDOM  SEARCH 

At  this  stage  it  seems  appropriate  to  apply  the 
mathematical  tools  developed  so  far,  and  to  discuss  two 
very  simple  and  well-known  types  of  searches:   one  which 
will  be  called  the  "systematic"  search,  and  another 

which  is  commonly  known  as  "random"  search. 

5 
In  the  systematic  search  the  observer  is  assumed 

to  move  in  such  a  way  that  no  part  of  the  area  is  searched 

more  than  once.   In  terms  of  the  model  this  means  that 

at  every  time  step  a  new  square  is  searched  with  probability 

one.   Hence,  the  search  process  is  specified  by 


f  1,  k  =  0, . . . ,n-l 

ass(k)  =  \ 
ss       [  0,  k  =  n. 


(8) 


From  equation  (8)  it  is  immediately  obvious  that  the 
A(m)'s  are  degenerate  random  variables  and  that 

m,  m  =  0  , .  .  .  ,n 

n,  m  =  n,n+l, ....  (9) 


EAgs(m)  =  J 


The  distribution  of  the  time-to-detection,  D   ,  is  obtained 

s  s 

from  equations  (1)  and  (9)  as 


0,  m  <_  0 

m      , 

-,  m  =  1,  .  .  .  ,n 


1,  m  >  n,  (10) 


^Koopman  [Ref.  3]  considers  this  search  as  the  first  case 
in  the  search  by  parallel  sweeps,  but  does  not  give  it  a  name 
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and  the  detection  rate,  d   (m) ,  from  equations  (2)  and 


(9)  as 


d   (m)  =   P[D   =  m|D   >  m-1] 
ss         ss     '  ss 


,  m  =  1, . . . ,n.  (11) 


n-m+1 

The  random  search  assumes  that  the  observer  is  equally 
likely  to  search  any  one  of  the  n  squares  at  every  time 
step,  no  matter  whether  he  has  searched  that  particular 
square  before  or  not.   Thus  the  probability  of  searching 
a  new  square  at  any  one  time  step  is  the  ratio  of  the 
number  of  squares  not  searched  so  far  and  the  total 
number  of  squares.   But  the  number  of  squares  not  searched 
so  far  is  just  the  total  number  of  squares  minus  the  num- 
ber of  squares  searched  so  far,  and  thus 

PCX   (m)  =  1|A   (m-1)  -  k]  -  Sdt  . 
L  rs       '  rs         J    n 

Therefore  the  sequence  which  specifies  the  search  process 

for  the  random  search  is 

a (k)  =  ~^  ,  k  =  0,...,n.  (12) 

Here  it  is  convenient  to  derive  the  detection  rate 

first.   From  equations  (6)  and  (12)  one  obtains 

d   (m)  e  P[D   =  m|D    >  m-1] 
rs         rs     '  rs      J 

Ea   [A   (m-1)] 
rsL  rs      J 


n-EArs(m-l) 

n-A   (m-1) 

E[ Eg ] 

n      J 


n-EArs(m-l) 
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-  E[n-EA   (m-1)] 
n       rs     J 

n-EArs(m-l) 
=  1,  m  =  1,2,...     .  (13) 

This  result  together  with  equation  (7)  can  be  used  to 

get  the  distribution  of  the  time-to-detection  as 

F   (m)  e  P[D    <  m] 
rs         rs  —  J 

m     , 
=  1-  n   [l-i] 
i-1   n 


■{ 


0,         m  <_  0 

l-C1^)   ,  m  =  1,2,...     .     (14) 


Of  course,  the  above  result  could  have  been  obtained 
more  easily  from  the  following  direct  argument :   at  any 
time  step  the  probability  of  not  detecting  the  target  is 

„   -I 

— - —  ;  hence  the  probability  of  not  detecting  the  target 

at  the  first  m  time  steps  is  ( — - — )  ,  and  thus  D   has 

n     '  rs 

the  geometric  distribution  with  success  probability 
p=  —  which  is  also  the  detection  rate. 

Clearly  the  systematic  search  is  optimal  for  the 
model,  but  it  is  equally  obvious  that  such  a  search 
"without  overlap"  can  be  performed  only  under  ideal 
conditions.   It  requires  no  less  than  that  the  observer 
have  perfect  knowledge  of  his  position,  perfect  recollec- 
tion of  his  path,  and  perfect  ability  to  move  wherever 
he  wants,  at  each  time  step  of  the  search.   The  random 
search,  on  the  other  hand,  can  be  regarded  as  minimal 
for  the  model,  in  that  an  observer  with  the  desire  to 
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detect  the  target  can  hardly  do  worse.   The  significance, 
then,  of  the  systematic  and  the  random  searches  is  not 
that  they  may  be  accurate  descriptions  of  real  world 
searches,  but  rather  the  fact  that  they  can  serve  as 
boundary  cases  for  certain  types  of  searches  which,  for 
lack  of  a  better  term,  will  be  called  "purposeful" 
searches. 
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6.   PURPOSEFUL  SEARCHES 

A  purposeful  search  can  be  defined  vaguely  as  a 
search  which  lies  somewhere  in  between  a  random  search 
and  a  systematic  search,  or,  in  other  words,  as  a  search 
which  is  at  least  as  "good"  as  a  random  search  (and  not 
"better"  than  a  systematic  search;  this  condition,  however, 
is  satisfied  by  all  possible  searches  in  the  context  of 
the  model).   To  make  the  definition  precise  some  criteria 
for  "good"  have  to  be  chosen.   The  choice,  of  course,  is 
arbitrary  as  long  as  the  criteria  are  stated  in  terms  of 
quantities  which  are  of  significance  in  a  search. 

One  such  quantity  is  the  time-to-detection;  a  "short" 
time-to-detection  is  "good",  and  hence  a  search  could 
be  called  purposeful  if  its  time-to-detection  is  smaller 
than  that  of  a  random  search  or,  alternatively,  if  its 
expected  time-to-detection   is  less  than  that  of  a  random 
search.   Another  important  quantity  is  the  detection 
rate.   Equations  (11)  and  (13)  show  that  the  systematic 
search  has  a  strictly  increasing  detection  rate,  whereas 
that  of  the  random  search  is  constant.   This  suggests 
as  a  condition  for  a  purposeful  search  that  its  detection 
rate  be  non-decreasing.   Intuitively  this  means  that  the 
more  time  the  observer  has  spent  searching  for  the  target 
without  finding  it,  the  more  likely  he  is  to  detect  it 
at  the  next  step. 

In  this  paper  a  search  will  be  called  a  purposeful 
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search  if  and  only  if 

(1)  its  time-to-detection  is  not  greater  than  the 
time-to-detection  in  the  random  search,  and 

(2)  its  detection  rate  is  non-decreasing. 
In  subsection  3.3  it  was  shown  that  a  sequence 

<a(k),  k  =  0,...,n>  uniquely  specifies  a  search  process; 
hence  the  conditions  above  must  be  equivalent  to  some 
conditions  on  the  a(k)  sequence. 

Sufficient  conditions  in  terms  of  the  a(k)  sequence 
for  a  search  to  be  purposeful  will  be  discussed  presently 
Before  that,  however,  a  brief  digression  to  the  theory 
of  stochastic  orderings  seems  appropriate. 
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7.   STOCHASTIC  ORDERINGS 

The  theory  of  stochastic  orderings  deals  with  order 
relations  between  random  variables;  some  results  of  this 
theory  are  presented  here  to  facilitate  the  understanding 
of  the  proofs  which  follow.   Theorems  are  not  stated  in 
their  most  general  forms  but  only  as  general  as  required, 
and  are  called  lemmas. 

A  random  variable  X  is  stochastically  greater  than  a 

st 
random  variable  Y,  denoted  by  X  >   Y,  if  and  only  if 

PCX  >  z]  1  P[Y  >  z] 

for  all  real  numbers  z . 

Lemma  1 .   Let  X  and  Y  be  non-negative  random  variables  such 
that  X  >St  Y.   Then  EX  >  EY. 


Proof 


00  00 

EX  -  /   PCX  >  z]  dz  >_  /   PCY  >  z]  dz  =  EY. 
0  0 


st         st 
Lemma  2 .   Let  X-,  ^_   Y,  ,  Xp  >_   Y?,  and  assume  that  (X-,,Y,) 

is  independent  of  (X2,Y2).   Then  (X1+X2)  >_st  (Y]L+Y2). 

Proof:   Let  F2  and  G,  be  the  distribution  functions  of 

X2  and  Y, ,  respectively.   Then 

00 

PC(X1+X2)  >  z]  =   /  PCXX  >  z-u]  dF2(u) 


>_     f    P[YX  >  z-u]  dF2(u) 
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=   /  P[X2  >  z-u]  dG1(u) 


>   /  P[Y2  >  z-u]  dG1(u) 


=  P[(Y1+Y2)  >  z]. 


Lemma  3 •   Let  X  and  Y  be  geometrically  distributed  random 

variables  with  success  probabilities  p,  and  p„  respec- 

st 
tively,  and  assume  p,  <_  p?.   Then  X  >_   Y. 

Proof: 

PCX  >  k]  =  (1-Pl)k  >  (l-p2)k  =  P[Y  >  k]. 

Lemma  4 .   Let  A(m)  be  the  number  of  events  in  the  area 
acquisition  process.   Then  for  every  m  =  0,1,..., 
A(m+1)  ^st  A(m). 
Proof:   From  equation  (3)  it  follows  that 

A(m+1)  =  A(m)  +  X(m+1). 
Hence 
P[A(m+l)  > k]  =  P[A(m)+X(m+l)  >  k] 

=  P[A(m)  >  k]  +  P[A(m)=k,X(m+l)=l] 

>  P[A(m)  >  k]. 
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8.   SUFFICIENT  CONDITIONS  FOR  PURPOSEFUL  SEARCHES 


The  conditions  a  purposeful  search,  as  defined  in 
section  6,  has  to  satisfy  can  be  expressed  formally  as 

(1)  D    <St  D 

v  '   ps  —    rs 

(2)  dps(m+l)  >_  dps(m),  m  =  1,2,...,6 

where  the  subscript  ps  denotes  a  purposeful  search.   The 
following  propositions  provide  tools  to  establish  conditions 
on  the  a(k)  sequence  of  a  purposeful  search. 


Proposition  1.   Let  {<A,(m),  m  =  0,1, ...>,Z}  and 
{<Ap(m),  m  =  0,1,... >,Z}  be  search  processes  specified  by 
the  sequences  <a-,  (k) ,  k  =  0,...,n>  and  <ap(k),  k  =  0,...,n>, 

respectively,  and  assume  for  k  =  0,...,n  that  a-,  (k)  <_  a«(k). 

st 
Then  D,  >_   Dp,  when  D,  ,  Dp  are  the  respective  times-to- 

detection. 

Proof:   Let  T1(k)  and  Tp(k),  k  =  l,...,n,  be  the  inter- 
event  times  of  the  processes.   By  definition  they  are 
geometrically  distributed  random  variables  with  success 
probabilities  a, (k-1)  and  ap(k-l),  respectively.   Since 
by  assumption  a-,(k)  <_  a?(k)  for  all  k  =  0,...,n,  thus 
Tl(k)  lSt  T2(k),  k  =  l,...,n,  by  Lemma  3-   The  inter- 
event  times  of  each  process  are  independent  by  definition; 

k  k 

hence  E   T.  (i)  >_s  l      T0(i),  k  =  l,...,n,  by  successive 

i=l  x        i=l  d  k 

applications  of  Lemma  2.   But  by  definition  1      T(i)  =  S(k), 

i  =  l 


In  the  case  of  the  systematic  search:  m=l,2, . . . ,n-l . 
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4-  i_ 

and  thus  the  k —  event  times  of  the  processes,  S, (k)  and 

st 
Sp(k),  are  stochastically  ordered  as  S,(k)  >_   S?(k), 

k  =  1,.  .  .  ,n,  and  hence  P[S1(k)  <_  m]  <_  P[S2(k)  <_  m],  by  defi- 
nition, for  all  real  numbers  m  and  integers  k  =  l,...,n. 

n  n 

From  this  it  follows  that   z   P[Sn  (k)  <_  m]  <_  i      P[S9(k)   m], 

k=l   1  k=l   d 

and  further,  using  equation  (1),  that  P[D,  <  m]  <_  P[Dp  <_  m] 

for  all  real  numbers  m.   Hence  by  definition  D,  >_   D?.  /~7 


Proposition  2.   Let  (<A(m),  m  =  0,1,. ..>,Z}  be  a  search 

process  specified  by  the  sequence  <a(k),  k  =  0,...,n>. 

Then  if  a(k)  is  linear  in  k  for  k  =  0,...,n-l  and 

a(n-l)  >_  — ,  the  corresponding  search  has  a  non-decreasing 

detection  rate. 

Proof:   From  the  assumption  of  linearity  it  follows  that 

for  k  =  0,...,n-l,  a(k)  =  a-bk  for  some  numbers  a,b. 

Since  for  every  search  process  a(0)  =  1,  hence  a  =  1; 

also  b  <  —  by  the  assumption  a(n-l)  >  — . 
—  n  J  *  —  n 

Let  b(x)  e  1-bx,  0  <_  x  <_  n-1,  0  <_  b  <_  — ;  then  for 

k  =  0,...,n-l,  6(k)  =  o(k).   Define  e(x)  =  ? ( X  \    , 

v  n— x ) 

0  <  x  <  n-1;  since  e'(x)  =  (1"nb^  >  0,  hence  e(x)  is 

"   "  (n-x)^  ~ 

non-decreasing  in  x. 

In  equation  (6)  it  was  shown  that  the  detection  rate  d(m) 

can  be  expressed  as 

d(m)  =  EaCA(m-l)] 
n-EA(m-l) 

for  all  integers  m  for  which  it  is  defined.   Since  in  this 
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expression  A(m-l)  can  only  take  on  the  values  0,...,n-l, 
it  can  be  re-written  as 

rffnO    EB[A(m-l)] 
d(m)  "  n-EA(m-l) 

or,  since  6  is  a  linear  function,  as 

_,  ,    B[EA(m-l)] 
d(m)  =  n-EA(m-l) 

which  is  the  same  as 

d(m)  =  6[EA(m-l)]. 

st 
In  Lemma  4  it  was  proved  that  A(m+1)  >_   A(m)  for  all  non- 
negative  integers  m;  thus  by  Lemma  1,  EA(m+l)  ^_  EA(m)  for 
m  =  0,1,...,  i.e.  EA(m)  is  increasing  in  m.   But  9  is  an 
increasing  function  of  its  argument,  and  hence  d(m)  is 
non-decreasing  in  m.   /~7 

From  Propositions  1  and  2  it  can  be  concluded  immediately 
that  if  the  sequence  <a(k),  k  =  0,...,n>  of  a  search  process 
is  such  that  for  k  =  0,...,n-l, 

(1)  o(k)  >  ars(k), 

(2)  a(k)  is  linear  in  k, 

then  the  corresponding  search  is  a  purposeful  search. 
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9.   FURTHER  DISCUSSION  OF  PURPOSEFUL  SEARCHES 

The  results  of  the  previous  section  seem  to  indicate 
that  any  search  with  non-decreasing  detection  rate  also 
has  a  time-to-detection  which  is  stochastically  less  than 
that  of  a  random  search.   That  this  is  indeed  the  case  can 
be  seen  from  the  following  argument . 

Consider  a  search  with  non-decreasing  detection  rate 
d(i),  i  =  1,2,...,  and  time-to-detection  D.   From  equation  (2) 
it  follows  that  for  any  search,  d(l)  =  — ;  hence  by  the 

assumption  of  non-decreasing  detection  rate,  d(i)>_— , 

m  i  m 

i  =  1,2,...   .   But  then  n   [l-d(i )  ]  <_  [1--]   for  positive 

i  =  l  n 

integers  m,  and  thus  by  equations  (7)  and  (14), 

Fn(m)  >  F   (m).   Hence  D  <st  D   . 
D    —  rs  —    rs 

Therefore,  clause  (1)  in  the  definition  of  a  purposeful 
search  is  redundant  and  can  be  omitted.   The  definition 
then  reads : 

A  purposeful  search  is  a  search  whose  detection 

rate  is  non-decreasing. 
Clause  (1)  could  now  be  added  as  a  theorem: 

The  time-to-detection  of  a  purposeful  search 

is  stochastically  less  than  the  time-to-detection 

of  a  random  search. 
The  sufficient  conditions  on  the  a(k)  sequence,  of  course, 
remain  unchanged. 
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10.   EXTENSIONS  OF  THE  SEARCH  PROCESS 

In  this  paper  it  was  shown  how  a  search  can  be  ex- 
pressed mathematically  as  a  stochastic  process  of  a  special 
type,  and  how  formulas  for  the  time-to-detection  and  the 
detection  rate  can  be  derived  from  this  process.   Because 
of  the  severe  limitations  on  the  model  under  consideration, 
this  effort  must  be  regarded  as  an  expository  example 
only.   To  obtain  results  which  are  of  theoretical  inter- 
est and  practical  use  it  is  necessary  to  extend  the 
search  process  to  more  general  models  which  take  into 
account  situations  where 

(1)  the  observer  uses  a  detection  device  which  is 
not  of  the  definite  range  law  type  and  which  has 

a  non-zero  false  alarm  probability, 

(2)  the  target  is  not  uniformly  distributed  and/or 
not  stationary, 

(3)  several  observers  search  for  a  target,  either 
independently  or  in  a  coordinated  action,  and 

(4)  the  observer (s)  and  the  target  move  in  continuous 
time  and  space. 

Some  of  these  extensions  seem  to  pose  no  difficulties. 
A  non-zero  probability  of  false  alarms,  for  instance, 
can  probably  be  dealt  with  by  stochastically  increasing 
the  time  between  successive  acquisition  events.   Other  ex- 
tensions will  certainly  require  a  great  amount  of  further 
research. 
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In  view  of  the  structural  similarities  between  a 
search  in  which  a  target  is  eventually  detected  and  a 
piece  of  equipment  which  eventually  fails,  and  subse- 
quently the  analogy  between  a  time-to-detection  and  a 
time-to-failure,  a  detection  rate  and  a  failure  rate, 
it  is  felt  that  this  research  could  benefit  by  using 
some  of  the  methods  and  tools  developed  in  reliability 
theory . 
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